Two-photon nonlinear spectroscopy of periodically trapped ultracold atoms in a cavity 
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We study the transmission spectra of a Bose Einstein condensate confined in an optical lattice 
interacting with two modes of a cavity via nonlinear two-photon transition. In particular we show 
that a nonlinear two-photon interaction between the superfluid (SF) phase and the Mott insulating 
(MI) phase of a Bose-Einstein condensate (BEC) and the cavity field show qualitatively different 
transmission spectra compared to the one-photon interaction. We found that when the BEC is in 
the Mott state, the usual normal mode splitting present in the one-photon transition is missing 
in the two-photon interaction. When the BEC is in the superfluid state, the transmission spectra 
shows the usual multiple lorentzian structure. However the separation between the lorentzians for 
the two-photon case is much larger than that for the one-photon case. This study could form the 
' basis for non-destructive high resolution Rydberg spectroscopy of ultracold atoms or two-photon 

spectroscopy of a gas of ultracold atomic hydrogen. 

; 

^ . PACS numbers: 32.80.Wr,42.50.-p,42.50pq 

INTRODUCTION 

^ I Cold atoms in optical lattices exhibit phenomena typical of solid state physics like the formation of energy bands, 
Josephson effects Bloch oscillations and strongly correlated phases. Many of these phenomena have been already the 
^ object of experimental investigations. For a recent review see Standard methods to observe quantum properties 

^ . of ultracold atoms are based on destructive matter- wave interference between atoms released from traps Q . Recently, 
' a new approach was proposed which is based on all optical measurements that conserve the number of atoms. It was 
shown that atomic quantum statistics can be mapped on transmission spectra of high-Q cavities, where atoms create 
a quantum refractive index. This was shown to be useful for studying phase transitions between Mott insulator and 
superfluid states since various phases show qualitatively distinct spectra Q. 

Experimental implementation of a combination of cold atoms and cavity QED (quantum electrodynamics) has 
made significant progress 0, Theoretically there have been some interesting work on the correlated atom- 

' O ' field dynamics in a cavity. It has been shown that the strong coupling of the condensed atoms to the cavity mode 
^ : changes the resonance frequency of the cavity Finite cavity response times lead to damping of the coupled 
atom-field excitations The driving field in the cavity can significantly enhance the locahzation and the cooling 
properties of the system[^ It has been shown that in a cavity the atomic back action on the field introduces 
atom-field entanglement which modifies the associated quantum phase transition 0, 11 1 • The light field and the atoms 



^ , become strongly entangled if the lattej: are in a superfluid state, in which case the photon statistics typically exhibits 

' complicated multimodal structures El2| . A coherent control over the superfluid properties of the BEC can also be 

. achieved with the cavity and pump [13| 

■ In this work, we show that a nonlinear two-photon interaction between the superfluid (SF) phase and the Mott 
insulating (MI) phase of a Bose-Einstein condensate (BEC) (confined in an optical lattice) and the cavity field show 

' qualitatively different transmission spectra compared to the one-photon interaction. Two-photon spectroscopy played 

I a very important role in the studies of BEC of atomic hydrogen [l^. Two-photon excitation of *^_R6 atoms to a 

■ Rydberg state was also achieved recently [1? 



THE EFFECTIVE TWO-PHOTON TRANSITION HAMILTONIAN 

The system we consider here is an ensemble of N two- level atoms with upper and lower states denoted by 1 1 > and 
|0 > respectively in an optical lattice with M sites formed by far off resonance standing wave laser beams inside a 
cavity. A region oik ^ M sites is coupled to two light modes as shown in Fig.l. Fig.l shows two cavities containing 
the two modes oi and 02 crossed by a one-dimensional optical lattice confining the BEC. In the two-photon process, 
an intermediate level \i > is involved, which is assumed to be coupled to |1 > and |0 > by dipole allowed transitions. 
The manybody Hamiltonian in the second quantized form is given by 




FIG. 1: Schematic diagram of the setup. The BEG atoms are periodically confined in an optical lattice and are made to interact 
with two laser modes ai and 02 which are confined in two intersecting cavities. Mode 02 is transmitted and measured by a 
detector. 
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In the field part of the Hamiltonian Hf^ ai are the annihilation operators of light modes with the frequencies toi, wave 
vectors k;, and mode functions ui{r), which can be pumped by coherent fields with amplitudes rji. In the atom part, 
Ha, ^l'(r) is the atomic matter-field operator, is the s-wave scattering length characterizing the direct interatomic 
interaction, and Hai is the atomic part of the single-particle Hamiltonian Hi. The detuning between the atomic 
transition frequency and any one of the two modes is nonzero. Under these circumstances, the intermediate state can 
be adiabatically eliminated and the effective Hamiltonian of the two-level atom can be written in the rotating-wave 
and dipole approximation as 

Hi=Hf + Hau (2a) 

Hal = TT ^ -TT'^^ ^ ihgola^ aiUi{r)a2U2{r) - H. c] (2b) 
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Here, p and r are the momentum and position operators of an atom of mass and resonance frequency uja-, c"*", cr~, 
and <Jz are the raising, lowering, and population difference operators, go is the atom-light coupling constant assumed 
to be same for both the modes. 

We will consider nonresonant interaction where the light-atom detunings A = wi -I- a;2 — Wa are much larger than 
the spontaneous emission rate and Rabi frequencies g^axai- Thus, in the Heisenberg equations obtained from the 
single-atom Hamiltonian H\ Q, can be set to —1 (approximation of linear dipoles). Moreover, the polarization 
a~ can be adiabatically eliminated and expressed via the fields a\ and a-z- An effective single-particle Hamiltonian 
that gives the corresponding Heisenberg equations for a\ and ai can be written as -fficff = Hf -I- Ha\ with 



Hai^^ + V^ci(r) + ^ala\aia2\u,{r)\'\u2{r)\\ (3) 



Here, we have also added a classical trapping potential of the lattice, T4i(r), corresponds to a strong classical 
standing wave. Interestingly, we find that unlike in \S\ we find that the Hamiltonian Hai does not contain terms like 
ul{r)u2{r) or ul{r)ui{r) which gives rise to an optical grating. 

To derive the generalized Bose-Hubbard Hamiltonian we expand the field operator ^'(r) in Eq. ([T]), using lo- 
calized Wannier functions corresponding to V^i(r) and keeping only the lowest vibrational state at each site: 
\I'(r) — ^fLibiw{r — Ti), where bi is the annihilation operator of an atom at the site i with the coordinate r^. 
Substituting this expansion in Eq. ^ with Hai ([3]), we get 

M 

^Y.^\Hh\h-il (4) 
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where 



Jf> = j dvw{v - r,) {^-^-^ + Ki(r)) w{v - r,), (5) 
J,^j = j drw{r - r,)|Mi(r)nu2(r)pu;(r - r^-), (6) 

U = ^TTash^/ma j dv\w[r)\^. (7) 

The BH Hamiltonian derived above is valid only for weak atom-field nonlinearity (3]. We assume that atomic 
tunneling is possible only to the nearest neighbor sites. Thus, coefficients ([5]) do not depend on the site indices 
{Jf\ = Jq' and Ji\±i ~ J'^^), while coefficients ([5]) are still index-dependent. The Hamiltonian ^ then reads 



H = Hf + J^^N + .r^B + ^4a|aia2 f ^ J^,^n. 
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\<ij> / i=l 

where < i,j > denotes the sum over neighboring pairs, fii = blbi is the atom number operator at the i-th site, and 
B = X]f=i ^l^i+i + H.c. While the total atom number determined by iV = X^i^li is conserved, the atom number at 
the illuminated sites, determined by Nk = X^iLi '"■ij not necessarily a conserved quantity. 

THE TRANSMISSION SPECTRA 

The Heisenberg equations for ai, 02 and bi can be obtained from the Hamiltonian ^ as 

ai = -I < -I- 2-^(X Ji,ihi + X '^Jj^^^j)a2«2 ) ai+rn (9) 

i=l <i,j> 

0,2 = -i {uJ2 + + X! '^».i^l^j)'^i"i r "2 + '72 (10) 

i=l <i,j> 
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bi = --^ ( Jq + —^ala{aia2JiA + Uhi ) bi 
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-■^ I ^ + —^al,a[aia2Ji,i+i ] bi+i 



\ ^ -^alalala2A^-^l j (11) 

Equations [MTD represent the coupled light-matter wave equations which determine completely the dynamics of the 
present system. In Eqns. [TOlandfTTl the second and third terms in the parentheses correspond to the phase shift of the 
light modes ai and 02 due to two-photon coherence. Note that the term that describes scattering of one mode into the 
other is absent. We consider a deep lattice formed by a strong classical potential Vci(r), so that the overlap between 
Wannier functions in Eqs. © and © is small. Thus, we can neglect the contribution of tunneling by putting J'^' = 
and Jij = for i ^ j. Under this approximation, the matter- wave dynamics is not essential for light scattering. In 
experiments, such situation can be realized because the time scale of light measurements can be much faster than 
the time scale of atomic tunneling. One of the well-known advantages of the optical lattices is their extremely high 
tunability. Thus, tuning the lattice potential, tunneling can be made very slow [l7| . In a deep lattice, the on-site 
coefficients Ji^i © can be approximated as Ji^i = |ui(ri)p|u2(ri)P neglecting details of the atomic localization. Using 
ai(t) = aiexp{—iuJit), we obtain the stationary solutions of Equations and 1101 as 

= TTl -2 (12) 

^Etil"i(r)PI"2(r)Pn.4«2 
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4^2 = — — ^"2 (13) 

Where Ap = cui — uj2- We now assume the mode ai to be in a the coherent state, which enables us to consider 

the quantity oi as a c number. The intensity of the mode ai is then large and undepleted as compared to the weak 

mode (02). As a result in Eqn. I12|, one could neglect the influence of the mode 02 on ai by ignorning the quantity 
2 

l'^'i(r)P|u2(r)pnia2a2 with respect to k. This basically means that we are ignorning the phase shift of the 
mode ai. This can be achieved by keeping the detuning A large and the probe amplitude 772 small. This yields 

2 

4^2 = — 2 (14) 

{Ap-^Etil"i(r)PI«2(r)Pn,| + 

Following Eq. (|14p allows to express a\a2 as a function /(ni, tt-m) of atomic occupation number operators 
and calculate their expectation values for prescribed atomic states 1 5') . 

For the Mott state {fii) mi = Qi atoms are well localized at the ith site with no number fluctuations. It is represented 
by a product of Fock states, i.e. I^*)}!!! = Y[i=i l^*)* = l^i' li^i)^ with expectation values 



if {hi, nM))Mi = f{qi, Qm), (15) 

For simplicity we consider equal average densities (ni)Mi = N/M = n {{Nk)mi — nK = Nk)- 
In SF state, each atom is delocalized over all sites leading to local number fluctuations. It is represented 
by superposition of Fock states corresponding to all possible distributions of N atoms at AI sites: |^)sf = 
J2q^ qf_, V^WM^/v^g7!^7gjv^T|gi^ qm)- Expectation values of light operators can be calculated from 

(/(ni,...,riM))gF = -r^v Y] ,^ ' J in, Qm) , (16) 

^ qil...qMl 

qi,---,qM 




^p|6 Ap/6 



FIG. 2: Photon number in the mode 02 for SF state (a) and MI state (b) due to one-photon transition (thin line) and two- 
photon transition (tliick line). The separation between the lorentzians in the SF state due to two-photon case is 115 and that 
due to one-photon case is 25. In the Mott state for the one-photon case, two-satellite contour is observed reflecting normal 
mode splitting of the two oscillators < ao,i > coupled through atoms. This splitting is absent in the two-photon case. All the 
frequencies are scaled with respect to S. Here 5 — go/ A. Also N — A4 = K = 15 and k = 0.35. The mode amplitudes are in 
scaled units 771 = 1 and 772 ~ 0.6. 




FIG. 3: Photon number in the mode 02 for SF state (a) and MI state (b) due to one-photon transition (thin line) and 
two-photon transition (thick line). The parameters are same as in Fig. 2 except k = 1.25. The comb like structure for the SF 
case (a) is replaced by a broad lorentzian for both the single-photon and two-photon case. For the Mott case (b), however, the 
peaks are broadened and the two-photon single peak shifts i.e, the dispersion is reduced. 



representing a sum of all possible "classical" terms. Thus, all these distributions contribute to scattering from a SF, 
which is obviously different from (/(ni, nM))Mi with only a single contributing term. We now apply this formalism 
to compare the two-photon spectra with the one photon spectra calculated in fs] Eqn. (5). For travelling waves, we 
take |'Ui(ri)p = |w2(ri)P = 1 and Ylf^i |wi(ri)p|M2(ri)pni = Ylf=i 

We present the transmission spectra for the one-photon case(thin line) and two photon (thick line) case in the SF 
state in Fig. 1(a) and that in the Mott state in Fig. 1(b) for the case where k — O.S^^/A. Clearly for both the cases, 
the transmission spectra is a sum of lorentzians with different dispersion shifts. A comb like structure is seen if each 
lorentzian is resolved. However the separation between the lorentzians for the two-photon case {US) is much larger 
than that for the one-photon case(2(5).Such a spectra can be reproduced experimentally by repeated measurements 
over a long time scale so that a superposition of different atomic distributions(obtained by tunneling) contibutes to the 
spectra. In the Mott state, the difference between the one-photon interaction(thin line) and the two-photon interaction 
(thick line) is striking. In the one-photon case, two-satellite contour is observed reflecting normal mode splitting of 
the two oscillators < ao,i > coupled through atoms. This splitting is absent in the two-photon case. A probable 



explanation for this observation can be traced back to the Hamiltonian Hai of Eqn. [21 In the absence of any optical 
grating the atomic grating is not formed and as a result diffraction of one mode into the other is absent (absence of 
Bragg scattering). In the one photon case, the two modes are bahaving like two independent oscillators exchanging 
energy via the atoms. In the two-photon transition, there is no energy exchange between the two modes but energy 
exchange occurs between the atoms and the two modes (taken together). The system then behaves effectively as a 
two- level atom interacting with a single mode via one-photon transition. The results for the case k > ffg/A is shown 
in Fig. 2. As found in Si] , the comb like structure for the SF case is replaced by a broad lorentzian for both the 
single-photon and two-photon case (Fig2a). For the Mott case, however, the peaks are broadened and the two-photon 
single peak shifts i.e, the dispersion is reduced. 

From the experimental point of view, such a two photon excitation has already been achieved in Rb atoms 
[lit ^^'^ atomic hydrogen BEC fl4|. Using the setup described in ^15], ®^i?6-atoms are magnetcally trapped in the 
55'i/2i^ = '2.,mF = 2 state and produce samples from thermal clouds to BECs by means of forced rf evaporation. 
After this preparation, the atoms are subject to a two-photon Rydberg excitation (using 780.246nm and 480. 6nm 
lasers) via the 5P3/2 state to the 43S'i/2 state with a length of the square light pulses between 170ns and 2^s. To 
reduce spontaneous photon scattering, the light is blue detuned by 2tt x ASiMHz from the 5P3/2, F = Z level. Thus 
only one photon per 100 atoms is scattered for the longest excitation time. However in such experiments the excitation 
is locally blocked by the van der Walls interaction between Rydberg atoms. 

CONCLUSION 

We have investigated off-resonant collective light scattering from a Bose Einstein condensate trapped in an optical 
lattice interacting with two cavity modes via nonlinear two-photon transition. Measuring the transmission spectra 
allows one to distinguish between a nonlinear two-photon interaction and a one-photon interaction. Depending on the 
state of the BEC, we found that the two-photon spectra is different from one- photon spectra. We found that when 
the BEC is in the Mott state, the usual normal mode splitting present in the one-photon transition is missing in the 
two-photon interaction. When the BEC is in the superfluid state, the transmission spectra shows the usual multiple 
lorentzian structure. However the separation between the lorentzians for the two-photon case is much larger than 
that for the one-photon case. This method could be useful for non-destructive high resolution Rydberg spectroscopy 
of ultracold atoms or two-photon spectroscopy of hydrogen BEC. 
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